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Abstract 

Compute-and-forward (CF) harnesses interference in a wireless network by allowing relays to compute combi¬ 
nations of source messages. The computed message combinations at relays are correlated, and so directly forwarding 
these combinations to a destination generally incurs information redundancy and spectrum inefficiency. To address 
this issue, we propose a novel relay strategy, termed compute-compress-and-forward (CCF). In CCF, source 
messages are encoded using nested lattice codes constructed on a chain of nested coding and shaping lattices. 
A key difference of CCF from CF is an extra compressing stage inserted in between the computing and forwarding 
stages of a relay, so as to reduce the forwarding information rate of the relay. The compressing stage at each relay 
consists of two operations: first to quantize the computed message combination on an appropriately chosen lattice 
(referred to as a quantization lattice), and then to take modulo on another lattice (referred to as a modulo lattice). 
We study the design of the quantization and modulo lattices and propose successive recovering algorithms to ensure 
the recoverability of source messages at destination. Based on that, we formulate a sum-rate maximization problem 
that is in general an NP-hard mixed integer program. A low-complexity algorithm is proposed to give a suboptimal 
solution. Numerical results are presented to demonstrate the superiority of CCF over the existing CF schemes. 

Index Terms 

Compute-compress-and-forward, compute-and-forward, physical-layer network coding, wireless relaying, nested 
lattice codes, modulo, quantization 


I. Introduction 

I NTERFERENCE has been long regarded as an adverse factor for wireless communications until the 
seminal work of compute-and-forward (CF) |[I|. The main idea of CF is to harness interference by 
allowing relays to compute linear combinations of source messages, without even the knowledge of 
any individual source messages. Since the advent of CF in |[T1, various CF-based schemes have been 
investigated in the literature. Much progress has been made towards the understanding of fundamental 
characterizations of wireless relay networks ll2]|- lfT^ . 

A CF scheme usually employs nested lattice codes [fT9ll - [l22ll . where the codebook of a nested lattice 
code is defined as the lattice points of a coding lattice confined within the fundamental Voronoi region of 
a nested shaping lattice If^ . [|23]|. In the original CF scheme [|T1|, a common shaping lattice is assumed 
for all source nodes, implying that every source node is forced to use a common power for transmission. 
This limits the potential of a CF scheme to exploit the asymmetry inherent in the nature of wireless 
communication channels. 

Recent work in 01, lISl presented modified CF schemes with asymmetrically constructed lattice codes, 
in which not only coding lattices but also shaping lattices are chosen from a chain of nested lattices. 
Precoding techniques based on channel state information (CSI) were proposed to enhance the performance 
of asymmetric CF schemes [|5||. However, both [JH and [l5]| were focused on asymmetry in the first hop 
of wireless relay networks. Not well understood is how the relays should process their received message 
combinations, or more specifically, how they should optimize the overall system when multiple hops of 
the network exhibit asymmetry. 
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Xiaojun Yuan is with School of Information Science and Technology, ShanghaiTech University, Shanghai, China. E-mail: 
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In a multi-hop relay network, the message combinations computed at relays are generally correlated, 
as they are generated from a common set of source messages. This implies that directly forwarding 
these combinations in general leads to information redundancy at destination. Meanwhile, the forwarding 
channel seen by each relay may vary significantly from each other due to the asymmetry of channel 
fading. It is thus desirable to reduce the forwarding rates of the relays with relatively bad channel quality. 
This inspires us to seek for more efficient relaying techniques for multi-hop relay networks. 

In this paper, we propose a novel relay strategy, termed compute-compress-and-forward (CCF). A key 
difference of CCF from CF, as manifested by their names, is an extra compressing stage inserted in 
between the computing and forwarding stages of each relay. The compressing stage at a relay consists 
of two operations: first to quantize the computed message combination on a lattice (referred to as a 
quantization lattice), and then to take modulo on another lattice (referred to as a modulo lattice). The 
design of the quantization and modulo operation should take into account the following two aspects. On 
one hand, it is desirable to choose a quantization lattice as coarse as possible and a modulo lattice as 
fine as possible, so as to minimize the forwarding rate at each relay. On the other hand, quantization 
and modulo operation in general suffer from information loss, and so the design of these compressing 
operations should be subject to the recoverability of source messages at destination. As such, there is a 
balance to strike in the design of the quantization and modulo lattices for compressing. 

To concretize the idea of CCF, we consider a two-hop relay network with multiple sources, multiple 
relays, and a single destination. The quantization and modulo lattices for compressing are respectively 
chosen as permutations of the coding and shaping lattices used for source coding. We present successive 
recovering algorithms to recover source messages at the destination. We then show that the above choice 
of compressing lattices is optimal in the sense of minimizing the forwarding sum rate under the constraint 
of the recoverability of source messages at the destination. Based on that, we formulate a sum-rate 
maximization problem that is generally an NP-hard mixed integer program. We propose a low-complexity 
suboptimal solution to this problem by utilizing the Lenstra-Lenstra-Lovsz (LLL) lattice basis reduction 
algorithm [I3]|, flM . Numerical results are presented to demonstrate the superiority of CCF over CF by 
exploiting the channel asymmetry of wireless networks. 

This paper is organized as follows. In Section |IIl we introduce the system model and some fundamentals 
of lattice and nested lattice codes. In Section Hill we describe asymmetric CF for the first hop of the 
considered network. In Section |IVl we describe the proposed CCF scheme involving quantization and 
modulo operation at relays. Section |Vl] is focused on the design of modulo operation, and Section |V] 
on quantization. The joint design of quantization and modulo operation is investigated in Section IVIII 
In Section IVIIIl we study the sum-rate maximization problem for the overall CCF scheme, and present 
numerical results to demonstrate the advantage of CCF over CF. Finally, the concluding remarks are 
presented in Section |IXl 


II. Preliminaries 


A. System Model 

Consider a relay network in which L source nodes transmit private messages to a common destination 
via M intermediate relay nodes. Each source node is equipped with a single antenna, and so is each 
relay node. Assume that there is no direct link between any source node and the destination. A two-hop 
relay protocol is employed. In the first hop, the source nodes transmit signals simultaneously to the relay 
nodes. In the second hop, the relay nodes transmit signals to the destination. The overall system model 
is illustrated in Fig. [H 

In the first hop, each source has a message w; G where is a finite field of size 7 and 7 is 

a prime number. Each source encodes w; as x; = / (w^) G and then transmit x; in the first-hop 

channel. The first-hop channel is a real Gaussian channel with additive white Gaussian noise (AWGN), 
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Fig. 1. A diagram representation of the system model. 


represented as 

L 

ym = '^hmlXl + Z^,m=l,--- ,M (1) 

1=1 

where ym G is the reeeived signal of the m-th relay, hmi ~ ^^(0,1) is the ehannel eoeffieient of 

the link from souree I to relay m, and is a Gaussian noise veetor drawn from AA(0,I„,). 

Denote by p/ = — ||x;|| the average power of souree 1. Then, the power eonstraint of souree I is given 
by 

Pi < Pi (2) 

where Pi is the power budget of source /. Further denote by H = [hmi] the first-hop channel matrix and 
by hm = [hmi, • • ■ , hmL]^ the channel vector to the m-th relay. 

In the second hop, each relay m communicates G to the destination. The second-hop channel 

is defined by the transfer probability density function p (y'|x'^, • • ■ , x'^), where y' is the received signal 
at the destination. The destination computes {w/} as an estimate of the original messages {w;}. Note that 
a detailed model of the second-hop channel is irrelevant to most discussions in this paper. Thus, we will 
only give an example of the second-hop channel later in Section IVTTTl 

For convenience of discussion, we henceforth assume rii = n 2 = n and L = M, i.e., the two hops have 
equal time duration and the number of the source nodes is equal to the number of relay nodes. 

We say that a rate tuple (ri, r 2 , • • • , r^) is achievable if 

Pr (w; = w;, / = 1, • • • , L) —)■ 0 as n ^ oo 

i.e., the destination can reliably recover the original messages {w;} by y' with a vanishing error probability 
as n —)■ cxD. This paper aims to analyze the performance of the network described above with CF-based 
relaying. 

B. Lattice and Nested Lattice Codes 

Nested lattice coding is a key technique used in CF-based relaying. To set the stage for further discussion, 
we introduce some basic properties of nested lattice codes. A lattice A C M” is a discrete group under 
the addition operation, and can be represented as 

A = {s = Gc : c e Z"} 
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where G G is a lattiee generating matrix [|20l . Let V denote the fundamental Voronoi region of A. 
Every x G ean be uniquely written as x = (x) + r, where r G V and Qx (x) is the quantization 

of X on A, i.e., the nearest lattiee point of x in A. Modulo-A operation [|^ is defined as 

X mod A = X — Qa (x) . (3) 

r llxlPda: 

The seeond moment per dimension is defined as (V) = „ vo^v)"^ ’ where Vol (V) is the volume of V. 
The normalized seeond moment of A is defined as G (A) = • We say that A is good for MSE 

quantization [|2T]| if 

lim G (A) = ^ (4) 

n^oo zne 

where e is the Euler’s number. 

A lattice Ai is nested in a lattice A 2 if Ai C A 2 . In this case, we say that Ai is coarser than A 2 , or A 2 
is finer than Ai. Eurther, if Ai C A 2 , then for x G M", 

[x mod Ai] mod A 2 = x mod A 2 . (5) 

A lattice codebook can be represented using a nested lattice pair (Ac, A^) with A^ C Ac, where A^ 
is referred to as a shaping lattice and Ac as a coding lattice. Denote the Voronoi regions of Ac and A^ 
respectively by Vc and Vs, and the corresponding volumes by Vc and Vg. The generated lattice codebook 
is 

C = Ac mod As = Ac n Vs. (6) 

The rate of this nested lattice code is given by 

i? =-log|C| =-log^. (7) 

n n Vc 

Moreover, we say that Ai, A 2 , • • • , Ak form a nested lattice chain if Ai 3 A 2 5 5 Ak lf25ll . Nested 

lattice codes with various rates can be constructed by appropriately selecting a pair of shaping and coding 
lattices from the chain, as detailed in Section |nll 

III. Eirst Hop: Asymmetric Compute-and-Eorward 

In this paper, we propose CCE for a two-hop relay channel, as illustrated in Eig. [2l This section is 
focused on the first hop which basically follows the asymmetric compute-and-forward (ACE) in Q, [[5]|. 
The ACE scheme involves asymmetric lattice coding with nested coding and shaping lattices, which 
improves performance by exploiting the knowledge of CSI. 



Fig. 2. The transceiver and relay operations for the two-hop relay network in Fig. [T] 
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A. Encoding at Sources 

We use nested lattiee codes to encode the messages of the sources. The lattices are generated following 


Construction A method in [(H, f[23 . Let k (•) be the mapping from the prime-sized finite field to the 
corresponding integers - ,7 — 1}, and (■) be the inverse mapping of k (•). Note that k (■) can 

be applied to a vector or matrix in an entry-wise manner. 

We first construct a chain of nested coding lattices following the Construction A method in [[Tl, [l22l . 
Let G G be a random matrix with i.i.d. elements uniformly drawn over F.^. Denote by Ga,i the 

first kA,i columns of G, where kA,i < fc is an integer. Define Ca,i = |G^^;b : b e F^'^’*|, and construct 

a lattice Finally, construct a coding lattice as Aa,i = BA^,/, where B G 

is a lattice generation matrix. In the construction, we require > /ca ,2 > ■ ■ ■ > kA,L, and so the 

constructed lattices are nested as 1 D Aa ,2 5 ^ Aa^l- Similarly, we construct a chain of nested 

shaping lattices A^^i 3 Ab ,2 5 ^ Ab,l, based on the same matrices G and B, with parameters 

kB,l > ^5,2 > ■ • • > kB,L- 

The lattice codebook Q for each source I is constructed as follows. We designate for the Z-th source 
a coding lattice Ac,; = Aa,-k^{i), and a shaping lattice A^,; = Ab,tts(i)^ where permutations tTc (•) and (•) 
are bijective mappings from {I,-- - , L} to {I,-- - , L}. In constructing lattice codebooks, we require that 
As,; is nested in Ac,;, and thus, kB,Trs{i) < kA,-K^(i)- Then the codebook of source / is C; = Ac,; fi V^,;. Note 
that the coding lattices {Ac,;} and shaping lattices {As,;} constructed above are good for both AWGN lH] 
and mean square error (MSB) quantization [1211 . 

We are now ready to describe the encoding function at each source. Let fcc,; = kA,-K^(i), and kg^i = k b, 7 ,^( 1 )■ 
The Z-th source draws a vector w; over F..y with length {kc^i — ks,i), and zero-pad w; to form a message 
as 



W; = [0^- • , 0, W; , 0, • • e F^^ . (8) 


The l-th source maps w/ to a lattice codeword in Ci as 



By following the proof of Lemma 5 in [[Tl, it can be shown that 0; (•) is a one-to-one mapping, which gives 
an isomorphism between the finite-field codebook and lattice codebook C;. Then, we construct 

the signal as 



where d; G is a random dithering signal uniformly distributed in the Voronoi region Vs,; of As,;. 

From Lemma 1 of lilTl . x; is uniformly distributed over Vs,;. Then, the average power of x; is given by 



( 11 ) 


where E (x) denotes the expectation of x, and Pi is the power budget of source Z in (|2l). 


B. Computing at Relays 

We now consider the relay operations. From ([T]) and (fTOl) . each relay m receives 

L 



( 12 ) 


i=i 


and computes a linear combination: 


L 



1=1 


(13) 
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where ami, I = 1, • • ■ ,L, are integer eoeffieients. To this end, the relay first multiplies y^, by and 
removes the dithering signals, yielding 


(J 


L 

m T E 

Z=1 

L 

^ ^ {flml (^Z T d;) T Omi^l) T 

Z=1 

L 

^ ^ (®mZ (tz Qhg^i (tz d/)) T T atmY‘m 

1=1 

L 

T ^ ^ 9ml^l T Olm'^m 
1=1 


(14) 


where step (a) follows by (fT^ and the definition of 9mi = ctmhmi — ami', step (6) follows from ([3]) and 
(fTOl) . Then the relay deeodes Sm by quantizing over a quantization lattiee Af^m, yielding 

Qhf „^ (®m) • (15) 

Following [HI, we ehoose Af^m as the finest lattiee in {Ac,;, for I with ami 7 ^ 0}, i.e. 

Af^m = fine {Ac,;, for I with ami 7 ^ 0 } (16) 

Note that Sm in (fT3l) is an integer linear eombination of ti, • • • , together with some residual dithering 
signals. This implies that not only the relays but also the destination is required to have the knowledge 
of d; for dither eancellation. 

We now determine the rate eonstraint to ensure the sueeess of eomputation at relay m. In (IT4l) . Zm — 
J2i=i ^mi^i + c^m'^m is the equivalent noise. An error in eomputing dm oeeurs when the equivalent noise 
Zm lies outside the fundamental Voronoi region of A^ This error probability goes to zero, i.e. 


provided 


lim PriSm^ <5m| = 0 

n—¥oo L J 


/ \ ^2 
14 ,; = Vol (Vc,z) > ( 27re max I 


(17) 

(18) 


where + Ya=i {ctmhmi - amif Pi is the power of z^, and = [a^i, • • • , amhf ■ 

Let P = diag (pi,p 2 , • • • ,Pl), and = diag • • • , y^). By (|7]), ([11]), and ([l8|), the rate 

of the Z-th souree is given by 


n = 


< 


1 , V 0 l(V.,;) 

n °®Vol(Vc,;) 




Pi 


max 

Kpl.dfjil ^0 


- log 
2 ^ 


Xjp=l iS^mhmp amp) Pj ^ 

( \ 

Pi 


max 


\m:a^l 70 


Q!^ + 


(oirrJ^m i 


(19) 


Note that the rate expression in ([T9l) reduees to the rate in Theorem 5 of [[D by letting Pi = P 2 = ■ ■ ■ = Pl- 
We will show that, instead of fixing pi = Pi, allowing pi < Pi leads to a eonsiderable performanee gain. 
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We now optimize {am} to obtain better eomputation rates. Denote 

iS^m) 

By letting (am) = 0, we obtain an MMSE eoeffieient as 


= 




1 + 


P^h. 


Substituting "m in (fT9l) . we obtain 


< 2 


Pi 


max ipm (aT) 


= - log"* 
2 ^ 




mm 


Pi 


\ 


P2a„ 


(h^Pa™)" 

— \r ~^—F 

1+ p^h™ 


A ~ 
= ri 


(20) 


( 21 ) 


(22a) 


(22b) 


where (x)"^ = max (x, 0). To summarize, a eomputation rate tuple (ri, r 2 , • • • , tl) is aehievable, i.e., (fTTI) 
is met, in the first hop if ri < fi, for I = 1, • • - L. Note that (l22l) reduees to the rate expression in |[T1 
Theorem 2] by letting Pi = ■■■ = = P. 


IV. Second Hop: Forwarding to Destination 

The preeeding seetion is foeused on the deeoding operation at relays. In what follows, we foeus on 
how to forward the decoded combinations dm to the destination. 


A. Compressing at Relays 

As illustrated in Fig. [2l after computing, each relay m compresses 5m and forwards G to the 
destination. at different relays are correlated, as they are combinations of the same set of source 

messages. Thus, the relay’s re-encoding problem is a distributed source coding problem. Forwarding 
directly at the relays in general leads to information redundancy at the destination. 

We propose the following two operations for the relays to compress |5m|- First quantize each 5m with 
lattice Kd^m, i-e. 

= QKd,rr. • (23) 

Then, take modulo of each ujm over a lattice Ae,m, i-O- 

Vm = <^m mod Ae,m- (24) 


The Wm obtained in (l24l) is a lattice codeword in the m-th relay’s equivalent codebook generated by 
the lattice pair (Ad,m,Ae,m)- Thus, with the above quantization and modulo operations, the forwarding 
rate of each relay m is reduced to 


TD A 

■^m 



V0l(Ve,m) y 

yo\iVd,m)J ’ 


(25) 


where the forwarding rate Rm is the rate of x(^. As illustrated in Fig. ^ Vm is then encoded as = 
£m (vm) and forwarded to the destination, where £m (■) is the re-encoding function of relay m. 
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B. Decoding at the Destination 

The decoding at the destination consists of two steps: (i) to compute v^, m = 1, • • ■ , L, from y', and 
(ii) to recover {w;} from {v^}. Without loss of generality, let TZ be the capacity region of the second-hop 
channel specified by p (y'|x'^, • • ■ , Then, in step (i), the destination can compute with a vanishing 
error probability, provided that the forwarding rate tuple (-Ri, • • • , Rl) satisfies 

(Ri, • • • , Rl) G R. ( 26 ) 

The remaining issue is to recover {w^} from {v^} in step (ii). We will discuss the design of the 
modulo lattices {Re,m} and the quantization lattices {A^m} to guarantee the recoverability of {w/} in the 
subsequent sections. 

C. Further Discussions 

The rest of this paper is mainly focused on the design of {Rd,m} and {Ke,m}- On one hand, it is desirable 
to choose {Arfm} as coarse as possible and {Re,m} as fine as possible, so as to reduce the forwarding 
rates at relays. On the other hand, {Arf^m} cannot be too coarse and {Re,m} cannot be too fine, so as to 
ensure the recovery of the source messages at destination. We will elaborate the design of {A^m} and 
{Ae,m} that ensures the recoverability of {w^} from y'. 

For ease of discussion, we henceforth assume no error in relay computation (i.e., 5m = Smi'^rn) and 
destination computation (i.e. the destination can perfectly recover {vm} with y'). Then, recovering {w/} 
from {y'} is equivalent to recovering {w/} from {v^}, where 

Vm Q^d.m (^rn) I^od Ag^m (2V) 

is the error-free version of v^- Since = (pi (w;) is an isomorphic mapping, recovering {w;} from {v^} 
is further equivalent to recovering {t/} from {v^}. 

V. Quantization at Relays 

In this section, we focus on the design of the quantization lattices {Ad,m} to ensure the recoverability 
of {t;} from {vm}. For convenience of discussion, we assume symmetric CF (SCF), i.e., all sources have 
the same power pi = p and thus the same shaping lattice A^ / = A^, V/. Then, the modulo lattices at relays 
can be trivially chosen as Ag^m = Ag,Vm. Then, in (1271) becomes 

Vm = i^rn) mod Ag 

= (ti-QA.(ti-di))^modA,. (28) 


A. Asymmetric Quantization Approach 

We now consider the design of the quantization lattices {A^ m}- We choose Ad,m-, m = 1, • • ■ , L, to be 
a permutation of the coding lattices Ag^/, Z = 1, • • • , L, i.e. 

Rd,m Rdim) 1 ^ 1) ' ' ' ) A (29) 

where (•) is a permutation function of {1, • • ■ , A}. 

The reason for the above choice of quantization lattices is explained as follows. From (|7]), (1251) . and 
(l29l) . the forwarding rate of the m-th relay is 

^rn = (30) 

As 7r“^ (vTd (■)) is a permutation, we obtain To ensure that the destination is able 

to recover all the source messages, the total forwarding rate can not be less than This implies 

that we can not choose {Arf^m} finer than (|29|) . 
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We say that (■) is feasible if the destination ean fully reeover {t/} from {v^}- In general, the 
quantization introduees information loss at relay m, and therefore (■) in (l29l) may be 

infeasible. We will show that, in a multi-relay system, sueh information loss at a relay does not neeessarily 
translate to information loss at the destination. 


B. Heuristic Discussions 

We first introduee the following faetorization of t/: 


= t; mod Kg 

Qaa ,2 (t«) + t; mod modAs 
^ 1,1 

Qaa.3 {Qaa. 2 i'^l))+QAA,2 (I/) mod AA,3+tl,l 


mod A, 


^ 4 - 

2 


-M=i 


mod A, 


(31) 


where e A^i^^ n Va,ii+i, for /r = 1 , • • • , L, (with Va,l+i = V^). Eaeh is a representation of t/ in 
the lattiee eodebook A^,^ n Va,^i+i- Therefore, = 0 if the shaping lattiee Aa,^+i is not eoarser than 
the coding lattice of t;, i.e. 


= 


t 

^=TTc(l) 


mod A^. 


(32) 


Consider Z = 1, • • • , L}, i.e., all the representations of in Aa,i n Va, 2 - Only one of them, i.e., 
non-zero, since the finest coding lattice Aa,i appears in {Ac,/} only once. 

We now consider recovering ^ from {v^}. From (l29l) . the only that contains ^ is 

Compute 


^- 1 ( 1 ) mod A A ,2 


QA^^-i(^)^^a^-i(i)«(t;-<5As(I/-d/))jmodA,modAA 

L 


1=1 

L 


mod A ^,2 




1=1 


mod A^ 2 


(^^-\i)^-\i)K-\i),i 


mod A^, 2 - 


(33) 


From |[T1, ^ is recoverable from (l3^ provided that n ^ j 7 ^ 0. Then, the contribution 

of 1 can be subtracted away from {vm}. 

We next consider recovering nnd t^-i( 2)2 which only appears in and v^-i( 2 )- By 

following similar steps in (l3^ . we obtain two linear equations of 2 t^-i( 2 ) 2 - Provrded that 

this linear system has a unique solution, we can recover 2 2 - Continue this process until 

all {t/,^} are recovered. Finally, {t/} are reconstructed using (l3^ . 
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C. Successive Recovering Algorithm 

We now present a sueeessive recovering algorithm by formalizing the heuristic discussions in the pre¬ 
ceding subsection. We first introduce some definitions. Denote the coefficient matrix A = [ai, • • • , G 
^LxL ];natrix A to the corresponding matrix Q over F..^ using k(-). Specifically, define Q = 

(A mod 7) G where the (f, j)-th element of Q, denoted by qij, is given by qij = {aij mod 7). 
Define the residual source set at the j-th iteration as 

(34) 

Define the residual relay set at the j-th iteration as 

r<''> = {/k,(/)<j,/G{i,---,L}}. (35) 

Define the residual coefficient matrix G F^^-^ as the submatrix of Q with the rows indexed by 
and the columns indexed by . In the above, the superscript represents the j-th iteration. Define 
the effective lattice codebook for the j-th iteration as Ce'^ generated by the lattice pair (A^ij, A^^+i)- 
Define a mapping from F^ to Ce'^ as 

(p,) A (Gb)] mod A^j+i (36) 

where 





Denote by (•) the inverse mapping of (■). 

We are now ready to present the Successive Recovering algorithm for asymmetric Quantization operation 
(referred to as the SRQ algorithm), as presented in Algorithm [T] below. 


Algorithm 1 (SRQ algorithm) 

Input: {Vm}m=l 5 5 

Output: b, / = 1, • • • , L 
1 : Initialization: A- for m = 1, • • • , L. 

2 : for j = 1 to L: 

3: Form a matrix by stacking {\m mod A^j+ij j ^ row by row. 

4 : Recovering: Compute •A ^ construct {w^j} by de-stacking W^) row by row, 

and compute <A (w^j), for I G 

5: Cancellation: If j < L, compute 



( ' '' V 

vb+i>^ 

m ' 

Vm ~ QAd,,n 1 tlml b,Iij 


yesO) r='^c{i) J 


mod A^ 


for m G 

6: end for 

'l2ii=Trc{j) 


mod As, for l 


I," 


,L. 


(38) 


We briefly explain Algorithm [T] as follows. In Line[3l rnaps a lattice point in Ce'’ to an integer 
vector in F^; {A^j} and are related by (|29|) : A^ ^+i = A^. In LineS ^ is the inverse of 

qO) in F^^A In Line[51 the contributions of G p < j} are cancelled from {v^}. 

A sufficient condition to ensure the success of Algorithm [U (i.e., t/ = t/, Z = 1, • • • , L) is presented 
below. 
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Lemma 1: Assume that = 1, • • • , L, are of full rank over F^. Then the output of Algorithm [T] 

satisfies / = 1, • • • , L. 

Proof: The proof is given in Appendix ■ 

The following lemma states the existenee of at least one feasible vr^ (•), i.e., all are of full rank 
such that we can recover source messages correctly by Lemma [B 

Lemma 2: For given tTc (•), if Q is of full rank over then there exists a mapping (•) satisfying 
(l29b such that every residual coefficient matrix is of full rank over F^, Vj. 

Proof: The proof is given in Appendix El ■ 

Note that feasible (•) is in general not unique. Thus, we need to search over all feasible vr^ (•) in 
performance optimization, as elaborated later in Section IVTTTl 


VI. Modulo Operation at Relays 


The previous section is devoted to the design of relay’s quantization operations when the sources use 
asymmetric coding lattices. In this section, we focus on the design of the modulo lattices {Agm} to 
reduce the forwarding rates. As analogous to the treatment in the preceding section, we trivially choose 
the quantization lattices as A^^m = Vm. Then, in (1271) becomes 


V 


m 


6m mod Ae,m 
■ L 

^ (f ^ Qa^i (fz ^z)) 

. 1=1 


modAe,m- 


(39) 


A. Asymmetric Modulo Operation 

Symmetric modulo operations have been previously used in CF fll], dUl to reduce the forwarding rates. 
In the symmetric modulo approach, each relay m takes modulo of 5m over the coarsest shaping lattice 
Ab,l, i.e., Ae,m = As,l, for m = 1, • • • , L. Then the forwarding rate at relay m is given by 


R 


m 



Vo1(Vb.l) A 

Vol(V/,™)J 


(40) 


where is the Voronoi region of Af^m defined in (fT^ . Such a rate of Rm may easily exceed the 
maximum of the source rates, which implies information redundancy. Therefore, the symmetric modulo 
approach is generally far from optimal. 

To further reduce the forwarding rates, we propose an asymmetric modulo approach to take modulo 
over different lattices at different relays. Specifically, we assume that each relay m takes modulo of 5m 
over Ae^m, with Ag^m, m = 1, • • • , L, being a permutation of the shaping lattices A^Z = 1, • • • , L, i.e.. 


Ae,m A^ jre(m)) m 1, • • • , L 

where Tg (•) is a permutation function of {1, • • ■ , L}. Then we reduce the forwarding rate to 


R 


m 



Vol 

yol{Vf,m) J 


(41) 


(42) 


For a random choice of TTg (•) in (|4TI) . the destination may be unable to recover {t;} from {v^}. We 
say that TTg (•) is feasible if the destination can correctly recover {tj upon receiving {vm}. 
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B. Heuristic Discussions 

To recover {t;} from {v^}, the main idea is to convert the asymmetric system (with {v^} in (l24l) 
defined using different shaping lattices) into a series of symmetric ones (each with a common shaping 
lattice). We start with the following observation on 


■ L 

■■ y^amljtl - - d;)) 

.1=1 


mod Ae,mniod A b,i 


■ L 

^ ami (t; - i - d/)) 

.1=1 


mod Ab,i 


^ ^ ami^i 


1=1 


mod As 1 


' L 

'^Umi {ti mod As^) 

. i 


mod A^ 1 


(43a) 

(43b) 

(43c) 

(43d) 

(43e) 


where As,i is the finest shaping lattice used in the system, (I43bl) follows from (l39l) . and (I43cl) follows 
from ([5]). Note that (t; mod As,i) in (I43el) is a lattice codeword in the fundamental Voronoi region of 
A^ i. Thus, (I43el) represents an effective symmetric system where {t/ mod A^ i} are treated as source 
messages defined over a common shaping lattice A^ i, and every relay m takes modulo of its received 
combination J2i ^mi (t; mod A^^i) over A^^i. Considering all the relays, this is a system of L equations 
with L unknowns. It can be shown that the system has a unique solution provided that the coefficient matrix 
A = [ai, • • • , is of full rank over [HI. However, the recovered (t; mod A^^i), Z = 1, • • ■ , L, 
are in general not equal to t;, except t; mod Ab,i = t; for I = (1). This exception is because 17 ^ 71 ( 1 ) 

is a lattice codeword in the Voronoi region of A^,!, and thus 177 - 1 ( 1 ) mod A^^i = 177 - 1 ( 1 )- To recover other 
codewords, we cancel the contribution of the correctly recovered t 77 ji(i) from the combinations {vm}, 
and discard v^-i(i) (as it is not useful in recovering other lattice codewords). In this way, we obtain a 
residual asymmetric system, where the (l)-th source and the (l)"th relay are deleted. Then, we 
can recover t^- 1 ^ 2 ) a similar way, and so forth. Finally, we can recover all {t/}. 


C. Successive Recovering Algorithm 

To make the above heuristic idea concrete, we introduce the following definitions. In the i-th iteration, 
the residual source set is defined as 

5«^{Z|7r,(Z)>^,/e{l,---,L}}; (44) 

the residual relay set is defined as 

7 "W 4 {/[TTe (m) > i, m G {1, • • ■ , L]] ; (45) 

the residual coefficient matrix Q G 1 + 1 ) x(i 1 + 1 ) submatrix of Q with the rows indexed by 5^*^ 
and the columns indexed by In the above, the superscript represents the Ath iteration. Define 
the effective lattice codebook for the i-th iteration as Ce , generated by the lattice pair (Aa,i, Ab^i). Define 
a mapping from F!^ to Ce'^ as 

fjb) A (Gb)] mod A^^j (46) 

where b G F!^, and the first kB,i bits of b are all zero. Denote by (■) the inverse mapping of (•). 

We are now ready to present the Successive Recovering algorithm for asymmetric Modulo operation 
(referred to as the SRM algorithm), as detailed below. 
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Algorithm 2 (SRM algorithm) 


Input: {'Vm}m=V {^c,«};=!’ 

Output: iij = 1, ■ ■ ■ , L 

1 : Initialization: for all m = 1, • • • , L. 

2 : for z = 1 to L: ^ 

3: Form by stacking mod j ,m G T*'*\ in a row-by-row manner. 


4: 


5: 


Recovering: Compute set as the corresponding row of and 

compute ^ ^ ^ 

d 


Cancellation: If ?' < L, then for m G denote —Qa , (t 

TTg TVs {1) ^ _s^7r_g ^ (i) \ 

calculate 


d*+i) 







— a -ir\ 

1 

) 

m rriTTs [ij 


J 

_ 


_ 


(i) - 

mod A.™. 


-1(0 


(47) 


6 : end for 


Algorithm [2] is briefly explained as follows. In Line [H we construct a symmetric system based on the 
mod-A^ i operation. In Line H we recover source codeword from the (^)-th source. In Line 

[51 we cancel the contribution of from |vm\m G T^*+d| to obtain 

Lemma 3: Assume that Q*'*\z = 1, • • • , L, are of full rank over F.y. Then the output of Algorithm [21 
satisfies t; = t;, / = 1, • • ■ , L. 

The proof of Lemma [3] is given in Appendix O Note that in Lemma [3l for given Q and {pi}, 
is a function of T*'*\ and thus a function of vTe (•). Therefore, Lemma [3] gives a sufficient condition for 
the feasibility of TTg (•). The following lemma ensures the existence of a feasible vTe (•) that makes all 
full-rank. Similarly to (•), feasible vTe (•) is in general not unique. 

Lemma 4: For given vr^ (•), if Q is of full rank over F.^, then there exists a mapping vTe (•) such that 
every residual coefficient matrix ig of full rank over F^. 

Proof: The proof is given in Appendix [Dl ■ 

VII. Asymmetric Quantization & Modulo Operation 

We are now ready to consider the joint design of quantization lattices {Ad,m} and modulo lattices 
{Ae,m} by combining the results obtained in the preceding two sections. 


A. Heuristic Discussions 

To start with, we now assume a general setting: the sources generally have different coding and shaping 
lattices; the relays perform both asymmetric quantization and modulo operation respectively defined by 
([29[) and ([4T]) . Then, in ([27[) becomes 


= (^m) mod Ae 


QAd,^\^arnl{tl-QA^ftl-di))jmodAe,r. 


(48) 


v;=i 


With permutations (•) and vTe (•) for optimization, the forwarding rate of the m-th relay is given by 

Hm. 


VOI (Vg.TTeM) 
n Vol (Va..,m) 


(49) 
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As a chain of shaping lattices are involved, we generalize (|3TI) as 




'L+l 


-^<=1 


mod 


(50) 


where e n Vyi,^+i, for /i = 1, • • ■ , L, (with Va,l+i = Vb,i), and ti^L+i e Ab,i n V*,;. 

Our goal is still to recover {t;} from {vm}- To this end, we first take modulo of Vm on As,i, yielding 


\ra mod Ab ,1 


Q^d,m (E ^ml (t; - Qa,,; (t; - d;)) I mod Ab,i 
Q^ d,m. j mod Ab, 

Q^d,m (t; mod A^^i) j mod A^.i, 


z=i 


(51a) 

(51b) 

(51c) 

(5 Id) 


where (I51cl) follows by noting A^ i 3 A^ ;. It is not difficult to see that (I51dl) is a special case of (1281) by 
letting As = Ab,i and d/ = 0 in (1281) . Thus, we use the SRQ algorithm to obtain t/ mod A^^, Z = 1, • • ■ , L. 
Afterwards, we cancel the contributions of {t/ mod A^ i} from v^- Then, the resulting system is defined 
based on a common coding lattice Ab,i, different shaping lattices {A^^/}, together with the codewords 
{Qasi (fz)}- Thus, we use the SRM algorithm to further recover {Qasi (fz)}- Finally, we reconstruct 
{tz} as 

ti = ti mod Ab,i + (tj), Z = 1, • • • , L. 


Note that the details of the above cancellation and recovery process can be found in Appendix |El 


B. Successive Recovering Algorithm 

We now formally present the Successive Recovering algorithm for asymmetric Modulo and Quantization 
(referred to as the SRMQ algorithm), as shown in Algorithm [3l In Line [2] of Algorithm [3l is a lattice 

codeword in the m-th relay’s equivalent codebook generated by the lattice pair (A^.i, Ae,m). 


Algorithm 3 (SRMQ algorithm) 


Input: {Ac,{^s,^}z=l’ 

Output: L, Z = 1, • • • , L 

1 : ^ SRQ({v^ mod , {Ac,;}^!^, , Ab,i). 

2 : ^ Vm- (Ez=i mod Ae,^, for m = 1, • • • , L. 

3: ^di- for Z = 1, • • • , L. 

4: 4- SRM({vr"}Lp (driztp {A5 ,i,---,A5,i}, {A,,,}f^^, 

{Ae,m}m=l- 

5: t; ^ + if°^, for I = !,■■■ ,L. 


Lemma 5: Assume that and j = 
Algorithm [3] satisfies L = t/, Z = 1, • • • , L. 
Proof: The proof is given in Appendix |El 


1, • • • , L, are of full rank over F^. Then the output of 
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C. Achievable Rates of the Overall Scheme 

We are now ready to present the aehievable rates of the proposed seheme. We have the following 
theorem. 

Theorem 1: For given H, A, and TZ, a transmission rate tuple (ri,r 2 , • • • ,rL) is achievable if there 
exist vTd (•) and vTe (•), such that the following conditions are met: 

1 ) power constraints in (fTTI) : pi < 

2) computation constraints in (I22bl) : ri < ri,\/l, 

3) forwarding rate constraints in (l26l) : (i?i, • • • , Rl) G 71, where Rm, m = 1, • • • , L, are given by 


Rm. 


r^-1 


TTc Pdim)) 


+7:log| 


P TTs ^(7re(m)) 
Pn~^{TTaim)) 



(52) 


4) recovery constraints: all and are of full rank over 

Proof: Condition [T] is the power constraint in (fTTI) . Condition [2] ensures a vanishing computing error 
at every relay; see (l22l) . Also, Rm in (15^ is obtained by noting (01), (|7]), (fTTI) . (|4T]) . (I29l) . and (1491) . Finally, 
Condition |4] ensures the recoverability of the source messages using the SRMQ algorithm as specified in 
Lemma [51 Therefore, (ri, r 2 , • • • , r/,) is achievable under Conditions fH to 01 ■ 

The following theorem ensures the existence of feasible (•) and vTe (•) that make all and 
full-rank. 

Theorem 2: For given tTc (•) and tt* (•), if Q is of full rank over F.^, then there exists mappings (■) 
and VTe (•) such that all residual coefficient matrices and are of full rank over F..^,. 

Proof: Lemma and Lemma 01 ensure the existence of feasible (■) and TTg (•) such that and 
Vi,j, are of full rank. ■ 


VIII. Sum-Rate Maximization 

In this section, we consider optimizing the achievable sum rate for the proposed CCF scheme. A 
centralized node is assumed to acquire all the knowledge of H, 71, and {Pi}. This centralized node 
informs each source I of lattice pair (Ac,;, and each relay m of quantization lattice Ad^m and modulo 
lattice Ac „i- 


A. Problem Formulation 

Based on Theorem fH we formulate the sum-rate maximization problem for the SRMQ algorithm as 
follows: 


max 

s.t. 


L 

1=1 

Pi < Ph'^i, 


( 


A < -log"^ 




mm 


Pi 


llp^aJb _ 

V " " i-pp^hJ^y 


,V/, 


(-Rl, • • • , Rl) G 7Z, 


Rm = r- 


1 , 


TTc ^(7rd(m))A2^^® 


Pit a ^(we(m)) 


\PtTc ^(TTa(m))^ 

rank = L — i + l,i = 1, ■ ■ ■ ,L, 

rank(Q<^>) = j, j = 1, • • • , L. 


(53a) 

(53b) 

(53c) 

(53d) 

(53e) 

(53f) 

(53g) 
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Note that if the SRM algorithm is used for message reeovery, then in (I53el) should be replaeed by 

= max (n + ^log , ( 54 ) 

y 2 y Pi J J 

Similarly, if the SRQ algorithm is used for reeovery, then (I53el) should be replaeed by (l30l) . 


B. Approximate Solution 

The problem in (l5^ is an NP-hard mixed integer program. Here we present a suboptimal solution as 
follows. For given {pi}, we apply the LLL algorithm [[3]|, [|2^ to determine the integer matrix A. For 
given A and {pi}, we seareh over all permutations (•) and all permutations vTe (•) to maximize (I53al) . 
This ean be done by noting that, given A, {pi}, (•)> and vie (•), the problem in (15^ is eonvex (as the 

eapaeity region TZ is always a eonvex set). Also not that the feasible set of (■) and vie (•) is not empty, 
as ensured by Theorem [2l The above diseussions are based on fixed {pi}. Finally, we need to maximize 
(I53al) by exhaustively searching over the quantized values of {pi}. 

We now give more details on determining A using the LLL algorithm. We determine A in a column- 
by-column manner from ai to a/,. For any step m, we first note that maximizing the right hand side of 
(153cl) is equivalent to minimizing pm (I22al) . We rewrite pm 

/ P^b P\ 

(II- ^ J Pa^ = a^D,^a^ (55) 

V l-f||Ph^|| / 

where 

Bm = P (il - ^ 

V i + ||Ph,„||V 

is symmetric and positive definite. Cholesky decomposition of gives = L^L^, and then 
is the Gram matrix for a lattice with generator matrix L^- Apply the LLL algorithm to to find the 
reduced matrix L(„. Compute A^ = Then we choose a^ as the row of A^ with the smallest 

norm that is linearly independent of a/, / = 1, • • • , m — 1. This ensures that the constructed integer matrix 
A is of full rank over F.^. 

The above algorithm involves an exhaustive search over {pi}, vrg (•) and vrrf(-). Thus, it is still com¬ 
putationally intensive, especially when the network size is relatively large. The study of more efficient 
algorithms for solving (15^ is out of the scope of this paper. 


C. Numerical Results 

To keep the computation complexity tractable, we set the second-hop channel be parallel channels as 

y'm = 9m^'m + ^'m 

where gm ^ ^ is the channel coefficient from relay m to the destination, g^, ~ AA(0,1), e M-xi 
is the forwarded signal at the m-th relay, 2 ||x(„||^ = PR^m, and G is i.i.d. Gaussian noise, 

~ AA (0, I„). So the rate region 7?. is a hypercube in M^. 

In simulation, the following settings are employed: L = 2; Pi = P,\/l; Abnite = 100; Pr,™, = 0.25P, Vm. 
Note that Abmte represent the number of discretized power levels in exhaustively searching each pi. 
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SNR (dB) 


Fig. 3. Performance comparison in a two-hop system. 

1) Comparison of Different Approaches in a Two-Hop System: We now compare the following schemes: 

• SCF: the original CF in [[T]|: 

• SCF-Q: SCF with asymmetric quantization (using the SRQ algorithm for recovering); 

• ACF: ACF with symmetric modulo and conventional quantization; 

• ACF-M: ACF with asymmetric modulo and conventional quantization (using the SRM algorithm); 

• ACF-MQ: ACF with asymmetric modulo and asymmetric quantization (using the SRMQ algorithm). 

The simulated sum-rate performances are compared in Fig. [3l We see that ACF outperforms SCF by 

about 2 dB. Also, SCF-Q outperforms SCF by about 2 dB, due to the use of asymmetric quantization; 
ACF-M outperforms ACF by about 2 dB, due to the use of asymmetric modulo operation; ACF-MQ 
outperforms ACF by about 3 dB, due to the use of both asymmetric modulo operation and asymmetric 
quantization; ACF-MQ outperforms SCF by about 5 dB, thanks to the use of asymmetric shaping lat¬ 
tices, asymmetric modulo operation, and asymmetric quantization. Fig. [3] demonstrates clearly the the 
performance advantage of the proposed CCF scheme over the conventional CF scheme. 

2) Performance Comparison With Different Network Sizes: We now provide the performance com¬ 
parison with different network sizes. Because of high computation complexity, we only simulate the 
performance of SCF and SCF-Q schemes. In Fig. HI we see that, for a 2 x 2 network (with two sources 
and two relays), the power gain of SCF-Q over SCF is about 2 dB at high SNR; for a 4 x 4 network, 
the corresponding power gain is about 6 dB. We see that the performance gain of SCF-Q increases 
significantly with the network size. 

It is highly desirable to understand how the performance of CCF improves with the network size. This 
requires the development of more efficient algorithms to solve (15^ . which is nevertheless an interesting 
topic for future research. 


IX. Conclusions 

In this paper, we proposed a novel relay strategy, named CCF, to exploit the asymmetry naturally 
inherent in wireless networks. Compared with conventional CF, CCF includes an extra compressing stage 
in between the computing and forwarding stages. We proposed to use quantization and modulo operation 
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Fig. 4. Performance comparison with different network sizes. 


for compressing the message eombinations eomputed at relays, whieh reduces the information redundancy 
in the messages forwarded by the relays and thereby improve the speetral effieieney of the network. 
Partieularly, we studied CCF design in a two-hop wireless relaying network involving multiple sourees, 
multiple relays, and a single destination. We derived aehievable rates of the seheme and formulate a sum- 
rate maximization problem for performanee optimization. Numerieal results were presented to demonstrate 
the signifieant performanee advantage of CCF over CF 


Appendix A 
Proof of Lemma [H 


We follow Algorithm[T]step by step to prove that the output of Algorithm[T]satisfies = t;, / = 1, • • ■ , L. 
Reeall that Algorithm [T] involves L iterations. In the j-th iteration, G is eomputed using 

|vm^, m e |, where Vm is set as for j = 1, and given in (l38l) for j > 1. 

We first prove that, in the j-th iteration. 



■ L L 

. 1=1 ^i=j 


mod 


(57) 


for m E . We prove this statement by induetion. In initialization, we have E For 

j = 1, the residual relay set only eontains one element, i.e., = {tt^^ (1)}, and = Aa,i is 
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the finest lattiee in {Ac,;}. Then 




(t; -d;)) 


1=1 

L 


mod Ao 


Ai); X] 


a-i 

1=1 iJ,=i 


mod A, 


(58a) 

(58b) 

(58e) 


where (I58bl) follows from (| 2 ^ and J2f=i (i)i ~ Q^s ~ d;)) G Aa,i; (I58el) follows from (ISTl) . 
Thus (157]) holds for j = 1. 

Now suppose that (l57l) holds for the j-th iteration. We show that tij = tij, I G in the j-th iteration, 
and that (1571) still holds for the (j + l)-th iteration. 

In Line [3l with < Vm \, we obtain 


v^^modA^j+i 


J=1 fJ,=j 
L 


1=1 


modA^j+i 
mod Aaj +1 


(59a) 

(59b) 






mod A^j+i 


(59e) 


for all m G where (I59al) holds since (1571) is true for the j-th iteration; (I59bl) follows from mod 
Aaj+i = 0, for /i > j, since G AA,ij.r\VA,ij.+i', (I59cl) follows from tij = 0 for / ^ by the definition 
of tij and Also by the definition of tij, we have tij G Ce'^ = Aaj D Vaj+i. In Line [3] of Algorithm 
[H maps t; j from Ce'^ to the finite field F^. Denote wij = (fzj)> for / = 1, • • • , A. From (l36l) 
and Lemma 6 in [[U, we obtain the following isomorphism: 


Then by (l59l) we have 


Qral^lJ = i’ 






mod Aaj+i 


X] = 'ip mod A^y+i) . 

ze5(‘) 


(60) 


(61) 


Stacking the two sides of (IbTl) for m G in a row-by-row manner gives a matrix equation over F.^: 

q6>w6 ) = U6> (62) 

where the rows in are specified by wfj, I G and TjO^ is the matrix obtained in (l38l) of Algorithm 
[TJ Recall that is of full rank by assumption. We obtain from (l62l) that 

wO) = (qO))”^ u6 > = W^'A (63) 

Thus, we can obtain wij, I G S^^'> from (l63l) . Then 

(w;j) = tij = iij, I G 56 ). 

Note that wfj is in general not the /-th row of W^A 


(64) 
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,(i+i> 


^ I ^ ^ Q 1 1 Q ^d,m O-ml ^l,/. 

1=1 V=1 // \leS(i'> fJ-=Trc(l) 

L L \ L / j 

QAd.mj ^ j ^fl mlQAs{^l ~ dO~QAd.m ( ^ ^ ^ml ^ ^ 

Z=1 /i=l / i=l 

^^d,m I ^ ^ ^ml ^ ^ I Q^d,7n O'ml U,fj, 

1=1 M=1 / \IgS^ 1'> /l=1Tc{l) 

L j \ L L / j 

Qf^d,m IE ^ml ®mZ E tz,/i Q^d,m E OmZ 


i^e^G) ii=wcii) 

mod Ao 


mod A. 


mod A, 


=1 n=l 
L L 

_z=l At=i+1 


z=i zi=i+i 
mod As 


Y^< 50 ) fl —' 7 rc(/) 


mod A, 


(66a) 

(66b) 

(66c) 

(66d) 

(66e) 

(66f) 


Then we eancel the eontribution of |X]|^= 7 rc(z) ^ ^ | from {v^, m G yielding (l66l) for 

m G where (I66bl) is from (1^ and (l28l) : (I66cl) and (I66el) utilizes the fact that, for xi G M”, X 2 G A, 


Qa (xi + X 2 ) = Xi + X 2 - (xi + X 2 ) mod A 

= Qa (xi) + X 2 


and (I66fl) follows by noting 


Q^d.m ( y~~! ^rnl f Z,/J j —QA^^rn I y~~! ^rnl tz,/. 


Z=1 


^l=l 


iZG5<J> lL=TTc{l) 


(65a) 

(65b) 


by the definition of S^^'> in (IMl) and = 0 for /x < 7ic{l). Eq. (l66l) shows that (1571) holds for the 
(j + l)-th iteration as well. Thus, (l57l) holds by induetion. 

Since Algorithm [H reeovers {tij, I G by (l64l) for j = 1, • • • , L, and tij = 0 for / ^ we see 
that all {t; j} are reeovered. Finally, Algorithm [T] reeovers {t;} by (l32l) . whieh eoneludes the proof. 
Furthermore, from (l64l) and t; ^ = 0 for I ^ we see that 


tz 


L 


L 


mod As = 

y~! 

_M=’rc(i) 


fl=TTc{l) 


mod As = ti 


whieh eoneludes the proof. 


Appendix B 
Proof of Femma[2] 

By assumption, = Q is of full rank. From (1^ and (1^ . we see that can be obtained by 

deleting one row (eorresponding to the (j)-th relay) and one eolumn (corresponding to the (j)-th 
souree) from 
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Suppose G > 1, is of full rank, and we want to find a G of full rank. 

Sinee tTc (•) is given, we delete the eorresponding eolumn of yielding a matrix of rank (j — 1). 
From linear algebra, there always exists at least one full-rank obtained by deleting one row of 

Choose sueh a (with rank (j — 1)) and set the eorresponding value of vr^ (•) aeeording to the 

index of the deleted row. By induetion, we obtain full-rank , and 


Appendix C 
Proof of Lemma [3] 

We follow Algorithmic] step by step to prove that the output of Algorithm [2] satisfies t/ = t;, Z = 1, • • • , L. 
Algorithm [2] involves L iterations. In the z-th iteration, is reeovered with |vm\m G where 

Vm is the lattice equation from the m-th relay at the z-th iteration. 

We first prove that, in the z-th iteration. 



ami (t/ - Qa,,; (t/ - di)) 

lesa) 


mod Ae,m 


(67) 


for m G and that exactly one of {t;} is restored in each iteration. We prove this statement by 

induction. The algorithm sets Vm^ = Vm,z7z G at the initialization stage. We immediately see that 
(1671) holds for z = 1. Now suppose that (1671) holds for the z-th iteration. We will show that 
recovered in the z-th iteration, and that (IFTI) also holds for the (z -f l)-th iteration. 

By the definition of and j = is the finest in {As /|Z G j, and also the finest in 

{Ae,m|’^ e T^)}. Then, with |vm we obtain 


v6) mod As,i 


^ ^ ami (fi Q^a,l 


mod Ab^i 


(t; mod AB,i) 


mod AB,i 


( 68 ) 


where mod A^ j is a codeword in the effective lattice codebook Ce'^ defined in Subsection IVI-CI That 
is, mod As j can be seen as from a common codebook Ce \ Note that 


since AB,^ = 

Define = ■^“6) (t^ mod Ab^i), for / = 1, • • ■ , A. By (1^ . (|4^ . and Lemma 6 in [HI, we have 

= z/;“6) (v6) mod AB,i) . (70) 

lesa) 

Form the matrix W^) by stacking w I G 56 ) , row by row. Similarly, form the matrix by stacking 
^Vm mod As^ij j ,m G We can write (iTOl) as 

Q(dw6) = U6), (71) 

which is a matrix equation over F!^. Then 

WW = (q6))-^u6) = W6) 


(72) 
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where is of full rank by assumption. From (1^ . we have w -i,.-, = . By setting w^_i as 

s \ / TTs 'TTs 

the eorresponding row of in Line S we reeover the message as and the lattiee 

eodeword as 



That is, the lattiee eodeword of the ^ (^)-th souree is eorreetly reeovered. 

Then the destination caneels the contribution of t 7 r 7 i(i) in Line [5] of Algorithmic 

to obtain By noting = 5^*^ \ (z)} and 7"^*+^) = \ {vr"^ (i)}, we have 


v(*+i) = 

m 


ami (t; - Qa,^i (ti - d;)) modAp 

l(zS(i+l) 

for m G which establishes (l67l) by induction. This concludes the proof of Theorem [3l 


(73) 


Appendix D 
Proof of Lemma [4] 

By assumption, Qd) = Q is of full rank. From (l44l) and (l45l) . we see that Qd+i) can be obtained by 
deleting one row (corresponding to the (j)-th relay) and one column (corresponding to the (j)-th 
source) from Qd). 

Suppose Qd) G > 1, is of full rank, and we want to find a Qd+d ^ Qf 

full rank. Since tt* (•) is given, we delete the corresponding column of Qd), yielding a matrix Qd) of rank 
(L — i)._From linear algebra, there always exists at least one full-rank Qd+i) obtained by deleting one 
row of Qd). Choose such a Qd+d (with rank {L — i)) and set the corresponding value of TTg (•) according 
to the index of the deleted row. By induction, we can obtain Qd), • • ■ , Q(^) that are all of full rank. 


Appendix E 
Proof of Lemma [5] 


In Section IVTIl we have shown that the SRQ algorithm can be applied on {v^ mod A^ i} to recover 
{t/ mod As i}. We next show that the SRM algorithm can be used to recover (Qa^ i (t;)}- From LineU 
of the SRMQ algorithm, we have 


Lquan 

^l 


L 


-M=l 


mod As,i 


■ L 

+ ti^L+l 

-M=l 

t; mod As,i 


mod Ab,i 


(74a) 

(74b) 

(74c) 


where step (I74bl) follows from ti^L+i G step (I74cl) from (l50l) . In Line IC we cancel from 

{v^}, yielding 


quan 


Vm Q^d,m ^ amli[ j mod A^ m 

■ L 

ami (t/ - tp" - Qa,^i (t; - d/)) mod A 

-i=i 

■ L 

y^ami (QAB.i(t;) -Qa^, (QAfl.i(ti) - (d;)) modA^, 


e,m 


. 1=1 


(75a) 

(75b) 

(75c) 

(75d) 
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where (I75el) follows from (1^ and (1^ . (I75dl) from the faet of = t/ mod i. Note that in (I75dl) 
is aetually a speeial ease of (1^ with t; replaced by Qab i (^0’ d/ replaced by = d; — for 
Z = 1, • • • , L. Therefore, we can apply the SRM algorithm to yielding outputs = Qas i (t;)}- 

Finally, we obtain 

t; = = ti mod Kb,I + Qab,i i^i) = 

for Z = 1, • • • , L, which follows from ([3]) and (1741) . 
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